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Angular bias momentum offers significant stability augmentation for hovering flight ve- 
hicles. The reliance of the vehicle on thrust vectoring for agility and disturbance rejection 
is greatly reduced with significant levels of stored angular momentum in the system. A me- 
thodical procedure for bias momentum sizing has been developed in previous studies. This 
current study provides groundwork for experimental validation of that method using an 
experimental vehicle called the Dual-Spin Test Device, a thrust-levitated platform. Using 
measured data the vehicle’s thrust vectoring units are modeled and a gust environment is 
designed and characterized. Control design is discussed. Preliminary experimental results 
of the vehicle constrained to three rotational degrees of freedom are compared to simula- 
tion for a case containing no bias momentum to validate the simulation. A simulation of a 
bias momentum dominant case is presented. 

1. Introduction 

Attitude stiffness is an attribute of a rigid body system’s rotational dynamics in which the body resists 
rotation in one or more axes, and tends to restore the disturbed system to its original attitude. This attribute 
can be passively achieved by introducing a constant angular momentum bias into the system so that the 
body resists rotation along the directions normal to that of the bias momentum. When the bias is provided 
by an internal flywheel, so that the body is not required to rotate about its bias axis, we have a “dual spin” 
system. 

Spacecraft designers have used dual spin for improving attitude performance for years, 1 2 but the technique 
is also very attractive for aeronautical vehicles. In [1], and other studies, it has been seen that introducing a 
sufficiently large angular momentum bias along the vehicle’s yaw axis results in a system that strongly resists 
and corrects deviations from trim in the pitch and roll axes. This passive attitude stiffness can substantially 
reduce the effort and bandwidth required of the vehicle’s fluidic control effectors, such as thrusters, in the 
case of hovercraft, or trailing edge surfaces for fixed-wing vehicles. 

This study focuses on validation of the methods discussed in reference [1]. This earlier paper provides 
expressions for sizing an angular momentum bias given a disturbance environment bandwidth and Power 
Spectral Density (PSD), and the desired Mean Square Response (MSR) of the vehicle Euler angles and 
body rates. The hardware, the Dual-Spin Test Device (DSTD) and related support equipment, was in 
a preliminary state in 2003 and was described in reference [3]. Now the hardware has matured and is 
ready for both experimental verification of bias momentum sizing expressions and experimental validation of 
bias momentum sizing methods. Unfortunately, due to some hardware failures, complete validation of bias 
momentum sizing is not offered in this publication. One test condition is offered and compared to simulation, 
however. Final results will appear in a publication at a later date. 

The following material will introduce dual spin dynamics and review bias momentum sizing as described in 
reference [1] and the issues relevant to its validation. Descriptions of the DSTD hardware and mathematical 
models and related system are offered. Control laws are synthesized for the system over a range of stored 
angular momentum levels; comparability is discussed. Finally, the groundwork for the validation test is 
documented including preliminary results for some test conditions. 
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2. Bias Momentum Sizing Overview 


A short summary of dual-spin dynamics and bias momentum sizing from reference [1] is presented for 
the sake of readability and to frame complications introduced by the DSTD hardware. The following first 
introduces the DSTD equations of motion with some simplifying assumptions. Then bias momentum sizing 
is summarized in general with some measures of system performance. Finally, a measure of the cumulative 
control activity is defined. 

The full non-linear equations of motion represented in scalar form are written in terms of the body axes 
which are defined in figure 1. These are simplified by fixing the body axes to the vehicle eg. The resulting 
equations (1,2,3) show coupling between all axes due to both inertial and gyric effects. 

Tl = IlUll + UJ 2 [h B + U)s(l3 — I 2 )] ( 1 ) 

t 2 = I 2 uj 2 + uji[h B + — Ii)] (2) 

T3 = /3W3 + h B + {I2 — I l)wi<x>2 ( 3 ) 

where I is the vehicle’s inertia, h B is the bias momentum, u>i is the angular velocity of the i th axis, and r, 
is the torque on the i th axis. 



Figure 1. Dual-Spin Test Device CAD model with axes 

Angular momentum products in each equation introduce nonlinearity to the system. Making a few as- 
sumptions yields more physical intuition and helpful simplification. First, we can assume the bias momentum 
to be sufficiently large such that the system is gyrically dominant, 1 i.e. 

\\h B \\»\\M\ (4) 

Furthermore, h B from Eq. 3 is assumed to be zero since the momentum wheel is regulated at a constant 
spin rate. The vehicle is assumed to have inertial symmetry about the roll and pitch axes eliminating the 
coupling in yaw. These assumptions yield a physically intuitive linear model of the dual spin system. 


1 ~ 

-U) 2 h B + Tl 

(5) 

1 2^2 ; 

« LOih B + r 2 

(6) 


I3U3 ~ T3 

(7) 


Yaw-axis-oriented bias momentum dominance causes cross coupling of the pitch and roll axis. Control 
torques applied to the pitch axis produce effects in the roll axis, and vice versa. So, the pitch and roll axes 
are cross controlled in a gyroscopically dominant system as outlined in reference [7]. 

The bias momentum sizing approach in [1] relates the mean square response of the vehicle in pitch and 
roll to the bandwidth and variance of a stationary ergodic torque disturbance assumed to come from the 
environment. For simplicity’s sake, 1 the disturbance is assumed to be independent in the roll and pitch axes, 
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with constant power spectral density (S Tl for roll, S T2 for pitch) over a common bandwidth BW T . This gives 
an approximate mean-square torque disturbance 

s(t t t)k2(S° i+ S° 2 )BW t (8) 


where e denotes expectation. This leads to the sizing equations 


where 




4 rTr } 

2cI(BW) 7 


arctan 


- 1 


arctan 


+ 1 


a BW t 

X„ = — : 



c system viscous damping coefficient 
viscous damping ratio 
/ vehical rotational inertia matrix 

A 0 undamped precession frequency 

BW t cutoff bandwidth of the disturbance torque 
h b bias momentum 


(9) 

(10) 


Plotting the mean square response for a particular disturbance environment yields some insight, 
damping is small and x 0 < 1 then Eq. 9 becomes approximately 




1 


h %( l ~ x l) 


e[r 7 


1 If 

( 11 ) 


The key insight from equation 11 is that the mean square response varies with the inverse square of bias 
momentum. In summary, one should choose bias momentum high enough so the disturbance bandwidth is 
much lower than the precession frequency and the mean square response is small. 

Bias momentum can reduce the required amount of high bandwidth active control of a hovering dual-spin 
system if the momentum bias is sufficiently large. Bias momentum cannot, however, mitigate persistent DC 
disturbances. Low bandwidth aerodynamic control is necessary to maintain trim. The Cumulative Control 
Activity (CCA) quantifies the cumulative amount of thrust vectoring used to maintain trim. The CCA is 
a sum of the squares of the propeller speeds (cc/ an J , and vane angles ( 0 van e 7 ) with respect to their trim 
values, Cjfam and d vanej , over time. The CCA can also be broken into the vane and propeller contributions 
CCA V ane CCi ^Fanspeeds’ 


CCA(t) 



& farii ) H - E ( @ vanej d vanej ) ] tit 
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C C. Ay anes (t) — I 'y ) i^vanej @ vanej) dt 
J 0 j=1 

ft 4 2 

C CAp anS p ee d s (t) = / ^ ) (^/art, &fani) dt 

do 


(12) 

(13) 

(14) 


3. DSTD Hardware 

The Dual-Spin Test Device (see figure 1) is a VTOL thrust levitating platform equipped with a shrouded 
propeller thrust vectoring system and a momentum wheel that adds a significant angular momentum bias 
to the vehicle. For the purposes of this paper testing is performed on the vehicle in a gimbaled configuration 
where the vehicle is attached to a bearing that allows 3 degrees of rotational freedom. The vehicle contains 
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Figure 2. Dual-Spin Test Device Hardware and CAD Model With Axes 


Parameter 

Symbol 

Value 

Vehicle Mass 

m 

18.15 kg 

Vehicle Inertias 

II, h; I 3 

0.8, 0.8, 1.0 ( kgm . 2 ) 

System damping 

c 

0.1 

Length from eg to fan centerline 

leg — ► fan 

0.4669 m 

shroud inner diameter 

D shroud 

0.3429 m 


Table 1. Values for DSTD Parameters 


several on-board and ground-based sensors for control and test measurements. Physical parameters including 
vehicle inertias, dimensions, component weights, and thrust vectoring data are included in table 1. 

The thrust vectoring system diverts the shrouded propeller thrust vector via a pair of control vanes 
oriented radially from the vehicle center. The pair may operate symmetrically or asymmetrically. Symmetric 
and asymmetric vane angle along with differential thrust between fan units yield sufficiently large controllable 
torques about the pitch, roll, and yaw axes. The propellers are proportionally slaved to the one across the 
diagonal of the platform as in 


^ fan\ T ^ fans ~ ^ fan 2 T ^ fan± 


(15) 


This keeps the magnitude of the total angular momentum vector nearly constant during maneuvers. The pro- 
peller rotation directions are counter-rotating for any two shrouds on an edge of the platform. The shrouded 
propeller static data, polynomial models, and control allocation scheme are included in the Appendix. 

The momentum wheel, shown in figure 3, provides an angular momentum bias proportional to its spin 
rate and rotational inertia. The momentum wheel is an aluminum spoke and rim design. A brushless DC 
motor drives the wheel within a safe operating range of approximately 0-4000 RPM. The wheels rotational 
inertia is 0.0411 kg-m which yields up to 17 Nms of angular momentum. This upper angular momentum 
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limit is sufficient to maintain a bias momentum dominant state where the vehicle dynamics are dominated 
by gyric effects. 1 



Figure 3. DSTD’s momentum wheel 


The DSTDs sensor suite measures the vehicle’s fan speeds, vane angles, attitude, angular rates, position, 
and accelerations. The suite is divided between on-board and ground based sensors. On-board sensors 
include a Microbotics MIDG inertial measurement unit. It measures accelerations and angular rates on all 
three body axes. Hall effect sensors mounted on the shrouded propeller motors sense a magnet mounted 
on the motor shaft to yield propeller RPM. Potentiometers sense the vane positions. All on-board sensors 
provide data to the on-board PC104 computer running a Simulink® program compiled through Matlab ®’’ s 
Real Time Workshop® . A ground based vision system measures the vehicle position and attitude. This 
system includes a camera and a PC. The camera transmits images of a dot pattern mounted on the vehicle 
to dedicated hardware in the PC which calculates and subsequently outputs position and attitude. These 
measurements are transmitted to the on-board computer or a ground based data system at approximately 
50 Hz. 

A JR3® load cell measures forces and moments on all three axes during non-flight static and dynamic 
thrust vectoring model parameter identification tests as well as some of the gimbaled bias momentum sizing 
tests. The load cell transmits data to dedicated hardware in a PC computer that is either recorded there or 
transmitted to the on-board computer system. 

A gimbal joint (see figure 4) provides three degrees of rotational freedom. The attachment point is 
adjustable so that the DSTD may be mounted at or below the eg; some counterweights are necessary. The 
gimbal consists of a shaft through two bearings, providing yawing freedom, and a universal joint, providing 
roll and pitch freedom. 

Four nozzles are situated vertically, one at the end of each shroud, (see figure 5). The nozzles are 
are connected by hoses to a bank of solenoid valves in a ground based system. These are supplied with 
compressed air. The valves can only be open or closed; there is no throttling. Forces are produced at the 
shroud tip when a solenoid opens. Since, the vehicle is gimbaled the forces are canceled at the rotation point. 
Only pure torques are applied by the system. 

4. Bias Momentum Sizing Validation Experiment Groundwork 

The section gives an overview of the groundwork for the bias momentum validation test. The DSTD 
propeller 2 motor burned up during testing. This prohibited the collection of data for any test conditions 
having stored angular momentum. An experiment was performed for the no-bias-momentum test condition. 
This experiment is compared to simulation. Good correlation between the experiment and simulation provide 
further validation of the simulation at the other test conditions as well, and thus to the concept of using 
bias momentum to augment vehicle stability. In place of measured data a simulation of the vehicle with 17.2 
Nms of bias momentum is presented and compared to the simulation at 0 Nms of bias momentum. 

The simulation results of the following sections are based on the DSTD mounted on its gimbal. This 
simulation contains non-linear equations of motion and the static and dynamic thrust vectoring system 
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Figure 4. Dual-Spin Test Device gimbal joint 


equation found in the appendix. The control law used for closed loop simulations is an LQR type and is 
described below. 

Trade studies between mean square response of the vehicle and the bias momentum levels show bias 
momentum improves performance in certain regions of the frequency ratio x 0 , (see equation 10). To experi- 
mentally verify these trades, bias momentum levels must be varied over a variety of disturbance bandwidths 
and the corresponding vehicle angular response must be measured. The frequency x Q begins the process of 
selecting numerical values. This ratio separates the trade space into three important regions. If x 0 > 1 the 
precession frequency of the system dominates the mean square response. However, if x 0 < 1 the precession 
frequency is not excited by the disturbance. So, the mean square response is attenuated in proportion to the 
bias momentum level. Since this is the region most advantageous for stabilization we chose the maximum 
amount of bias momentum our momentum wheel can produce, approximately 17 Nms. The BW T is then 
limited to ensure x 0 < 1. 

A bench test of the disturbance system identifies the maximum variance of the disturbance torque at the 
max BW t . This yields the spectral densities ( S ° x , S ° 2 ) of the torque disturbance about the roll and pitch 
axes. The gust disturbance system is driven by a random binary signal of known bandwidth. Each thruster 
receives it own unique random signal. The supply pressure varies in proportion to the number of jets on at 
any time. This results in a torque about the pitch and roll axes having a constant power spectral density up 
to 1.5 Hz, the cutoff bandwidth. The PSD of this signal can be seen in figure 7. The air jets create torque 
time histories, as measured by the JR3 load cell, on the vehicle that may be seen in figure 6. The torque 
produced, seen in figure 7, has a nearly constant power until approximately 1.5 Hz. A structural resonance 
of the vehicle exists at 10 Hz. This can be seen in the PSD of the torque signal in figure 8. The vehicle will 
need to be structurally reinforced before testing with bias momentum to avoid excitation of the precession 
frequency from structural vibrations. 

The jet system may be run at different pressures and bandwidths. The result is a repeatable and scalable 
simulated gust environment of known power spectral density, bandwidth, and variance valid for static testing 
and hover conditions. 
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Figure 5. Gust simulation jet 


T~Gust 

MSR(Nm 2 ) 

BW T 

T x 

0.3141 

« 1.5 Hz 

Ty 

0.2335 

« 1.5 Hz 

T z 

0.0590 

« 1.5 Hz 


Table 2. Disturbance Torque Sizing Parameters 




100 110 



Figure 6. Torque time history of gust disturbance 
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Frequency (Hz) 


Figure 7. Power Spectral Density of air jet signal 



Figure 8. 


Power Spectral Density of gust disturbance torque 
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A multi-variable Linear Quadratic Regulator (LQR) control posed as a tracking problem is used to return 
the vehicles attitude to a trimmed state or track an attitude command. The controller is synthesized using 
nine states including three angular velocities, three Euler angles and the integrated error between the Euler 
angles and the trim value for all three axes. Control is synthesized using the same LQR weights for each 
bias momentum level to ensure fair comparison. The control effort is defined as the CCA found in equation 
12. The cost for the LQR problem is in terms of the attitude errors {Pi), the vehicle body rates (07), and 
the integral error of the attitude angles (7*) is 

J(t) = f ({Pi ~ Pi ) 2 + (wj — Wj) 2 + (7 i — Ti) 2 ^ dt, (16) 

Jo i= 1 V 7 

where a)j, pi, and Pi denote trim values for each parameter. 

The linear DSTD model (5,6,7) is used for control design. Physical quantities may be found in the 
appendix. The trim condition is straight and level, meaning that the attitude angles are all set to zero. The 
lifting force produced is less than that necessary to lift the platform which reduces wear on the hardware 
and provides a more relaxed testing environment. 

4.1. Experiment and Simulation Results 

At this point only one experimental test condition was performed. The test is run with no bias momentum 
in the system. The following figures (9,10,3,11) compare the simulated test with the measured data. The 
same torque disturbance is applied to both the simulation and the DSTD. Close correlation between the 
simulation and the experiment validates the simulation. 

The first two plots, in figures 9 and 10, show the Euler angles and rates for both the simulation and 
DSTD experiment. The simulation predicts the DSTD’s Mean Square Response (MSR) of the pitch and roll 
attitude angles to within 5.6% and the vehicle body rates to within 1.9% , (see tables 3 and 4). The yaw 
angle does not correlate as well between the simulation and experiment. The control allocation algorithm 
assumes that the yaw command for the vane angles may be superimposed over the asymmetric commands. 
The larger, yet still small, error here is probably a symptom of this assumption when the vanes are not 
operating in their linear region. This assumption may need to be revisited if large angle maneuvers are 
attempted. 

The power spectral density of the simulated and experimental cases attitude angles in figure 11 have a 
relatively flat power until 0.2 Hz where they drop off slowly till 1.5 Hz. They drop off rapidly at this point 
since the disturbance at this frequency rolls off rapidly as well, as seen figure 8. The Cumulative Control 
Activity of the vanes, seen in figure 12, shows that the simulation and experiment agree to within 7.3% 
at the worst point. The CCApanspeeds agrees poorly being off by a factor of 3.5 believed to be caused by 
excessive sensor noise in the measured data and some fan model errors in the simulation. This plot has three 
separate pairs of curves. The largest curve is the total CCA found in equation 12. The smaller curves are 
the contributions of the vanes and the propeller speeds to the total CCA. 


Euler Angles 

M S Rexp {Deg 2 ) 

MSRsim(Deg 2 ) 

Roll 

1.865 

1.660 

Pitch 

1.317 

1.710 

Yaw 

0.557 

0.032 

Veh. Body Rates 

MSR Exp (%£)* 

MSRsim 

Roll 

13.028 

13.727 

Pitch 

13.189 

12.977 

Yaw 

1.096 

0.245 


Table 3. Mean Square Response of Euler angles and vehicle rates for experimental and simulated data 
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Figure 9. Time history of Euler angles for simulation and experiment 
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Figure 10. Vehicle attitude rates for simulation and experiment 
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PSD of Vehicle Attitude h = 0 Nms 



Frequency (Hz) 


Figure 11. PSD of Euler Angles for simulation and experiment 



Time (sec) 


Figure 12. Cumulative Control Activity for simulation and experiment 
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4.2. Simulation Results for Large Bias Momentum Level 

Due to significant hardware failures on the DSTD the results for cases containing bias momentum cannot 
be presented until repairs are made. Instead, the advantages of bias momentum dominated systems are 
reiterated here with the current simulation. 

The following plots show the previous simulation of the DSTD under the disturbance environment in 
figure 6 with hg = 0 Nms. The maximum level of Hb = VTNms is also included. The PSD seen in figure 15 
shows the spectral power for the pitch and roll of the hs = 1 7 Nms case is far smaller than the hs = 0 Nms 
case over the low bandwidth region. The high bias momentum simulation maintains trim with no increase 
in control activity. The MSR, (table 4), of the attitude and rates shows an order of magnitude improvement 
over the no bias momentum case. The CCA in figure shows no increase between the two cases. 
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Figure 13. Euler attitude angles for simulated experiment between bias momentum levels 


Euler Angles 

M SRhB=i7Nms(Deg 2 ) 

M S RhB=oNms (Deg 2 ) 

Roll 

0.223 

1.629 

Pitch 

0.222 

1.677 

Yaw 

0.029 

0.031 

Veil. Body Rates 

MSR hB =nNms 

MSR hB =0Nms (%*)* 

Roll 

3.479 

13.463 

Pitch 

3.725 

13.727 

Yaw 

0.2398 

0.239 


Table 4. Mean Square Response of Euler angles and vehicle body rates for no bias momentum and large bias 
momentum 
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Figure 14. Vehicle body rates for simulated experiment between bias momentum levels 
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Figure 15. PSD for simulated experiment between bias momentum levels 
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Figure 16. CCA for simulated experiment between bias momentum levels 
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5. Conclusions 


Limited validation of bias momentum sizing for dual spin systems was presented. An experiment was 
conducted that applied a known disturbance to the Dual-Spin Test Device. This disturbance has known and 
reproducible cutoff bandwidth and spectral power. The vehicle’s Mean Square Response and Cumulative 
Control Activity were measured and compared to simulation. A PSD of the vehicle’s attitude response showed 
attenuation by the controller and the bandwidth cutoff of the disturbance environment. The simulation’s 
PSD, Mean Square Response, and Cumulative Control Activity agree closely to experimental values in the 
pitch and roll axes. 

A simulation was run of the Dual-Spin Test Device under a bias momentum dominant condition and 
showed the value of bias momentum for reducing the vehicle’s inherent reaction to gusts. Confidence in 
the simulation at this test condition is strengthened considering the agreement between the simulation and 
experiment at the no bias momentum test condition. A future experiment of the bias momentum dominant 
test conditions is suggested. 
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Appendix 

A. DSTD Thrust Vectoring Model 

The thrust vectoring model relates fan speeds and vane angles to net propulsive forces and torques on 
the vehicle. Static bench test data was used to generate an experimentally verified polynomial model of the 
system. The data presented here is tailored for producing pure torques for attitude control. 

A. 0.1. Symmetric Vane Configuration Static Data 

Symmetric vane configurations direct the vane trailing edges of a shroud in the same direction. These are 
used to create yaw torques. A bench test is run to identify the ability of the shrouds to create torques about 
the vehicle yaw axis. Steady state data is taken oven the operating range of the shrouds. Figure 17 below 
shows the measured data and a polynomial model fit to this data. 



Figure 17. Yaw torque as a function of vane angle and fan speed 


A. 0.2. Asymmetric Vane Configuration Static Data 

Asymmetric vane configurations occlude the flow through the shrouds by directing the vane trailing edges 
either toward or away from each other on a shroud unit. Occluding the flow on one shroud produces a thrust 
differential between it and the opposing shroud. This allows for pitch and roll torque creation. In figure 18 
the change in thrust on a shroud is shown. This data is a polynomial fit to bench test data plotted in 19. 
The polynomial fit in 18 is the data found in figure 19 with the value of the zero vane angle data subtracted 
from the other data points at every propeller speed. In this way the change of thrust from asymmetric vane 
angles is known instead of the thrust value itself. In this model 9 > 0 indicates vane trailing edges are 
rotated inward by the amount of theta, likewise 9 < 0 indicates vanes rotating away at angle 9. 
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Figure 18. Change in thrust using asymmetric vane configurations 
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Figure 19. Change in thrust using asymmetric vane configurations 
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A. 0.3. Fan Thrust and Differential Thrusting 

Opposing shrouded propeller units may also create pitch and roll torques about the vehicle eg by differential 
thrusting. A bench test was performed over a range of fan RPM to identify a nominal thrust model for a 
single propeller unit. Figure [20] shows the measured data and a second order polynomial model. 



Figure 20. Shrouded propeller thrust over nominal range of RPM values 


A.0.4. Control Allocation Configuration 

Combinations of symmetric and asymmetric vane angles, and differential propeller thrusting provide con- 
trollable torques about all three axes at two different band widths. The control vanes can operate with 
significantly higher bandwidth than the propeller motors. High frequency torque commands are sent to the 
inverse vane models via a high pass filter. The remaining low frequency content is sent to the differential 
thrusting inverse models. The bandwidth of the propellers and vanes are determined via dynamic tests the 
nominal operating conditions. The vanes are limited to a bandwidth of 10 Hz at the nominal thrust condition 
where the propellers run at 4500 RPM. Both the vane and motor dynamics are assumed to be of the first 
order. The transfer function found from measured data and used in simulation to follow. 


TF, 
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